(19) 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(43) Date of publication: 

03.01.2001 Bulletin 2001/01 



(ID EP 1 065 799 A1 

EUROPEAN PATENT APPLICATION 

(51) lntci7: H04B 1/707 



(21) Application number: 99305245.5 

(22) Date of filing: 02.07.1999 



(84) Designated Contracting States: 

AT BE CH CY DE DK ES Fl FR GB GR IE IT LI LU 
MC NL PT SE 

Designated Extension States: 
AL LT LV NIK RO SI 

(71) Applicant: LUCENT TECHNOLOGIES INC. 
Murray Hill, New Jersey 07974-0636 (US) 



(72) inventor: Kariml, Hamid R. 

Swindon, Wiltshire SN5 5BJ (GB) 

(74) Representative: Williams, David John et al 
Lucent Technologies UK Limited, 
5 Mornlngton Road 
Woodford Green, Essex IG8 0TU (GB) 



(54) Multi-user detection for wideband CDMA 



(57) A method tor detecting user symbols in a multi- 
rate multi-user environment, of a WCDMA communica- 
tion system, the method comprising altering the order in 



which information symbols are estimated so as to group 
user symbols in a block-banded structure. Specific al- 
gorithms for achieving this are disclosed. 



O) 
CD 

m 

CD 

o 



1 K t 



m 

in 



m 



A 

±J. 



JR. 



Q. 
LU 



Printed by Jouve. 75001 PARIS (PR) 



EP 1 065 799 A1 



Description 

Field of the Invention 

s [0001] This invention relates to communication systems. In particular, it relates to mobile communications in a mutti- 
rate environment using wideband code division multiple access (WCDMA) 

[0002] In the multi-rate environment of the WCDMA uplink, the sampled signal block received at the base^station 
detector can be written as a vector r = Ax + n where x is the vector of bits (or symbols) transmitted by all K mobile 
users in the cell, n is the noise vector and A is a sparse matrix which is a function of the user spreading codes and the 
10 radio propagation channel. The objective of multi-user detection (MUD) is to estimate the vector x of ail user symbols, 
given the received vector r and the matrix A. 

[0003] The computational complexity of MUD can be prohibitive for a large number of users and is furthermore 
dependent on the structure of matrix A. 

15 Description of the Prior Art 

[0004] A communication system such as a mobile telephone communication system, may typically have several 
users using the system simultaneously. The mobile stations may be synchronised or asynchronous. Figure 1 illustrates 
the uplink transmissions of a number K of synchronised mobile stations. In the absence of inter-symbol interference 

20 (or when the delay spread is much smaller than the symbol period) it can be seen from the figure the nth symbol of 
any user interferes with the mth symbol of another user if and only if m - n. The result of this is that joint detection of 
the user symbolsx 1 ... x K can be achieved at the base station receiver by considering one single symbol epoch at a time. 
[0005] In an asynchronous uplink scenario, however, the base station receiver needs to process signals transmitted 
by any number K uncoordinated mobile stations. At a symbol level, the maximum asynchronous delay x max , between 

25 the transmission of any two mobiles is less than a common-symbol epoch (C^ chip intervals). A control symbol epoch 
is chosen since this is common among all users. Figure 2 shows such a system and it is seen that the nth symbol of 
a user interferes not only with the nth, but also with (n± 1) th symbols of all other users. Accordingly, joint detection 
of the user symbols can only be performed optimally by considering a signal block of an appropriate length. The block 
needs to be many symbols long. 

30 [0006] A natural MUD block length for WCDMA transmission is the time-slot period. There are two good reasons for 
this. Firstly, the channel impulse response can be assumed to be fixed over a time-slot interval (625 u.s) resulting in 
simplifications and reductions in the complexity of the multi-user detector. Secondly, longer block lengths (eg a frame 
period of 1 0 ms) can result in unacceptable latencies. This is especially the case for the detection of some of the control 
bits which may be required within a time slot period, such as power control or rate information. 

35 [0007] Consider K asynchronous mobile stations, each transmitting time slots with relative delays of xW<ri®< ... ^tW 
"samples" with respect to an arbitrary user labelled user 1 

( X 0> = o by definition). As can be seen from Figure 2 the MUD block needs to be ti(N c Q c +W-1 )+ t<K> samples long in 
order to account for Inter-Symbol Interference (ISI) and asynchronous delays. ISI represents the 'spill-over 1 between 
consecutive symbols transmitted by a user due to delays in the propagation channel. Note that the block-length is time- 
40 varying and depends on the maximum delay xflQflriQo-l , and is therefore upper bounded by 
MUD Block Length ^([N C +1]Q C +W-1)-1 samples 

Let (m^eC^NcQ^w- 1 )** represent the MUD signal block at the output of the base station's antenna-m receive filter 
when sampled at t\ times the chip rate. This may then be modelled as 



< m) r = ~ (m Ux + (m) n 

so where 

(K) 

< m >A€C^< Nc ^ c+w - l)+T I** 1 

55 

is the multi-user convolution matrix at antenna - m, x = [x< 1 > T ... rf ( > 1 I T eR N is the vector of N = Z k N( k > (data and control) 
symbols transmitted by 

all users within one time-slot and additive noise at antenna-m is given by 
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(K) 

<m>neC n < N cQc +w * l >* T 

5 The nature of the 

(*>Ax 

10 product can be further clarified by an expansion of equation(1 ) which is shown graphically in Figure 3 for K asynchronous 
users . 

[0008] In Figure 3, <m)AMeCfo(NcQc+w-i)x(i<j j s the convolution matrix for user k at antenna-m as defined and 



x<k> eR N 

is the vector of all NW symbols transmitted by user k over a time-slot period. 

[0009] Consider now a two-user uplink scenario. Suppose that user 1 accesses a 2-code voice service transmitting 
20 N/ 1 > = N c = 10 control symbols, x/ 1 > = x c < 1 >, with a spreading of <V 1 > = Q c = 256 and N 2 < 1 > = N d < 1 > = 20 data symbols, 
*2 (1) = *d (1) « witn a spreading of 

Q 2 0) =*Q d ( 1 ) = 128. User 2 accesses a 3-code data service transmitting N,< 2 ) = N c = 10 control symbols, x-, (2) = x c (2) , 
with a spreading of Q/ 2 ) = Q c = 256, N 2 < 2 > = N^ 2 ) = 20 data symbols, X2 (2) = x d1 < 2 >, with a spreading of C^^) = Q d1 ( 2 ) 
= 128 and N 3 < 2 ) = N^ 2 ) = 40 data symbols, x 3 < 2 > = x^ 2 *. with a spreading of Q 3 < 2 > = Q^ 2 ) = 64. 
2$ [001 0] Furthermore, user 2 has a maximum asynchronous delay of = r|Q c -1 with 
respect to user 1. The product ^Ax is graphically illustrated in Figure 4. In the Figure, 

the first two groups of N c ( 1 > and N d 0) columns relate to the first user; N c < 1 > shows the control symbols and N d < 1 ) shows 
data symbols. The following three groups of N c < 2 ), N d1 < 2 > and N^ 2 ) columns relate to the second user and show, 
respectively, control symbols and the two sets of data symbols for user 2. The length of each vertical line represents 
30 Q+W-1 which is equivalent to the length Q of the spreading sequence when filtered by a channel of length W chips. 
The chip oversampling factor in this example is i\. 

[0011] The objective of a multi-user detection algorithm is to "reverse" the influence of the convolution matrix, MA, 
on the transmitted information. Accordingly, MUD can be computationally very complex, particularly for a large number 
of users and is dependant upon the structure of matrix < m >A. 
35 [001 2] The present invention arose in an attempt to reduce computational complexity for multi-rate multi-user detec- 
tion, in WCDMA (wideband code division multiple access) systems. 

Summary of the invention 

40 [001 3] According to the present invention there is provided a method for detecting user symbols in a multi-rate multi- 
user environment, of a WCDMA communication system, the method comprising altering the order in which information 
symbols are estimated so as to group user symbols in a block-banded structure. 

[0014] Typically, each user transmits control symbols and at least one group of data symbols, and wherein a convo- 
lution matrix < m )A is formed which is block-banded, each block including a number of columns relating to the number 
45 of groups of data or symbols transmitted by the users, and grouped according to the respective control signals and in 
an order related to the order at which symbols arrive at a receiver. 

[0015] If control symbols are not considered in the multi-user detection problem, then the data symbols with the 
longest duration transmitted by each user effectively take up the role of control symbols. Other grouping methods may 
be used. 

50 

Brief Description of the Drawings 

[0016] Embodiments of the invention will now be described, by way of example only, with reference to the accom- 
panying drawings, in which: 

55 

Figure 1 shows uplink transmissions, in a WCDMA system, of K synchronised mobile stations; 

Figure 2 shows uplink transmission of K asynchronous mobile stations; 

Figure 3 shows a convolution matrix for K users at the mth antenna (antenna-m); 



3 



EP 1 065 799 A1 

Figure 4 shows the convolution matrix for two users; 

Figure 5 shows a column-permutated multi-user convolution matrix in the form of a block banded structure; and 
Figure 6 shows a column-permutated convolution matrix for two users. 

s [0017] It is described above that the signal, < m >r, at the output of the antenna-m receive filter can always be repre- 
sented as a linear function of a multi-user convolution matrix, < m >A, the 
vector of user symbols,*, and additive noise < m >n. 

[001 8] The objective of a multi-user detection algorithm is to eliminate the effects of multiple-access and inter-symbol 
interference by "reversing 1 the influence of the convolution matrix, < m >A, on the transmitted information. As such , efficient 

to implementation of MUD relies heavily on the exploitation of the structure and sparse nature of < m )A. 

However, as was seen earlier, while (™>A is always sparse, this cannot be readily exploited since the detailed structure 
of < m )A is highly dependent on the number of users active in the cell, their relative delays and mix of data rates. 
This issue is addressed in embodiments of the invention by considering certain permutations of < m >A which render its 
structure more predictable and hence computationally friendly. 

is [001 9] Consider the multi-user convolution matrix < m >A of Figure 3, corresponding to a MUD block-length of one time- 
slot. The matrix has an irregular structure which is a function of the asynchronous delays between the mobile trans- 
missions, as reflected in the offsets between sub-matrices ("OAW where k=1 ... K. This structure is a direct result of the 
fact that the elements of vector x (unknown information symbols) are grouped together according to the mobile user 
from which they originate. A benefit of this ordering scheme is that the symbol estimates at the output of a MUD 

20 algorithm could be readily assigned to the different users. 

[0020] The inventor has realised that it is possible to obtain a more predictable structure for < m )A by changing the 
order in which the information symbols are estimated. 

This would be equivalent to a scrambling of the elements of vecbr x and hence the columns of matrix ( m >A as described 
by 

25 

(m)Ax = (m)AVWx = < m >Ax, 



30 where V is a unitary permutation matrix and Vtysl. Following such a permutation, the model of Eq. (1 ) can be written as 

(2) (m) r= (m) Ax + (m) n 

35 [0021] Inspection of the structure of ( m )A reveals that a suitable permutation is one that groups the user symbols 
according to their respective control symbols and in a fashion which reflects the order in which the symbols arrive at 
the receiver. Such a permutation imposes a block-banded structure upon matrix MA. 

[0022] Note that for any multi-code scenario or mix of data rates, it is possible to derive a column-permuted multi- 
user convolution matrix ( m )A which is block-banded. Furthermore, when representing a block-length of one time- 
40 slot, < m >A can be completely defined via sub-matrices MB, _ (m >B Nc each containing K c columns and (a maximum of) 
T|(2Q C +W-1)-1 rows, as shown in Figure 5. 
[0023] In the Figure, 

45 K M*> Qc 

k=l i=l Qi (k > 

[0024] The column permutation process may be described as follows. For every control-symbol epoch, user k syn- 
so chronously transmits a total of £W = z, Q^QjM (i=1 .. . MM) symbols (one control and data). As a first step, shuffle 
the columns of ( m )AWVk=1 ... K, such that those columns corresponding to the above symbols are placed adjacent to 
one another. This ensures that each of the MaW sub-matrices are block-banded, consisting of N c blocks (one block 
for each control symbol in a time-slot) each of £M columns. As a second step, shuffle the generated blocks across the 
boundaries of the (™)A( k >Vk=1 ... K so that the blocks corresponding to control symbols which arrive at the base station 
ss at an earlier time are placed in the left-most columns. This ensures that the total convolution matrix MA is block- 
banded, consisting of N c blocks each of 1^ £M = K c columns. 

[0025] The validity of the above concept can be illustrated by considering a specific example. 

[0026] An arbitrary scenario will be considered where K=2 users are active on the uplink. User 1 transmits with 
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spreading factors of 0,0) = 256 and Q 2 < 1 > = 128, while user 2 transmits with spreading factors of Q^ 2 ) = 256, Q 2 < 2 > = 
128 and Q 3 < 2 ) = 64, resulting in the convolution structure depicted in Figure (4). The structure of the permuted 
matrix < m >A, as illustrated by Figure 5, is graphically illustrated in Figure 6, which shows the structure of column-per- 
muted convolution matrix ( m )A=< m >AV for a K=2 user scenario and chip oversampling factor of r\. Note that the block- 

5 length represented corresponds to a WCDMA time-slot containing N c =10 control symbols. Note that while the value 
of Kc depends on the number of users and their mix of data rates, the overall structure of < m >A does not. Each block 
comprises a first control signal , from user 1 , two data symbols and D 2 from user 1 , a second control symbol C 2 , 
from user 2, two data symbols D 3 and D 4 from the Nd^ 2 ) data columns of user 2 in Figure (4), and four data symbols 
D 5 , D 6 , D 7 , D 8 from the Nd 2 ^ data columns of user 2 in Figure (4). 

10 [0027] For the general case of K users, the sub-matrix ( m >B may be calculated as follows: 

MB=0[ti(2Q c +W.1)-1,Kc] 

15 

io=l 

for k=l...K 

20 for i=L..M^) 

Qi<k) 

nw=Q c /Qi< k > 

25 

epochs} 

(m)Bi(k)=0[ti(r i C)Q i (k)+W.l),ri( k >] 

code of k lh user} 

30 

if code is "Inphase" then oc=ai< k >=0 
35 if code is "Quadrature" then a=ai< k) =l 

for n=l...n*) {«' h symbol} 

40 d=fGi»> {[Csc< k >] n <n-.)Q ,W >«:nnQ ^Ci*)}*™^) 

[(m)Bj(k)] M ( n .l)^< k > + l: n (nQ (k W.|).n=d 

45 end 

(m)B((m)x(k)+l : <™H*>+Ti(ri< k >Qi< k W-l),io : i 0 + rV»-l) = <"»Bi( k > 

50 

io = io + n< k >; 
end 

55 end 

where for convenience, the number of rows in < m >B is assumed to be at the maximum value of t](2Q c +W-1 )-1 implying 



{k lh user} 
{i th code} 

{Spreading factor} 

{No of symbols over Q c chip 

{sub-matrix of B for i lh 
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the presence of a user with a maximum asynchronous delay of = t]Q c -1 . 

[0028] Note that the permutation implicit in the above algorithm implies that the elements o 

of vector x are scrambled versions of the unpermuted symbol vector x = [x < 1 > T ... x^ T l T - 

[0029] Symbols XjM transmitted by the i* code of user k may then be extracted from the permuted symbol vector x 
5 as follows: 



70 



15 



20 



25 



30 



35 



40 



45 



50 



{Deal with symbols in the nth 
{Start at column h of block B} 



Zt=l 

for k = 1:K 

for n = 1 : N c 

copy of block B} 

fori=l:M< k > 
Qi< k) 

r i ^> = Q c /Q i (k) 

Xi*>( [n-1] n< k > + 1 : nTj«) = x( [n-l]Kc+h : [n-1] Kc+/atn«.l) 
h = k + 
end 

end 

l :x = 0 {h counts the symbols of user k in 

matrix B} 

for i = 1:Mq* k > 



ri« = Qc/Qi« 

end 

Ii = h + fe 



{Update starting point for next 



user} 
end 

[0030] The above are examples of suitable algorithms. Other forms of block-banded structu re than those specifically 
shown may be used. 

[0031] It is possible that control symbols are not considered in the multi-user detection problem; they may be detected 
by some other means. In such cases, data symbols of chosen characteristics (e.g. those of longest duration) take up 
the role of control symbols in the above schemes and algorithms. The invention extends to these methods. 
[0032] In this specification, the following notation and nomenclature is used: 
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10 



MATRIX NOTATION 

[0033] 

x 
X 

NOMENCLATURE 



Vector. 
Matrix. 

Transpose and Hermitian Transpose. 

Element of Xon the /*» row and 7 th column. 

Sub-matrix of Xfrom rows r^ to r% and columns rr^ to rr^. 



75 



20 



25 



[0034] 



A 
A 

(m)g 



Convolution matrix for I th code of A* user at m lh antenna. 
Convolution matrix for A* user at m th antenna. 
Un-permuted multi-user convolution matrix at m* 1 antenna. 
Un-permuted multi-user antenna array convolution matrix. 
Permuted multi-user convolution matrix at antenna. 
Permuted multi-user antenna array convolution matrix. 
6 £>X1Qc+wa)-\ x Kc sub-matrix of the multi-user convolution matrix {m) A for 171 th 
antenna {m) A. 



30 


K 


■ it 1 


Number of users (maximum = 52). 




K c 






3S 




e R l 


Number of codes associated with the A* user. 


N?> 




Number of symbols transmitted via I th code of A* user in a 








time-slot. 








Number of symbols transmitted by the A th user in a 


40 






time-slot. 




N 


e* 1 


Number of symbols transmitted by all users in a time-slot. 




1 




Oversampling factor (nominally =2). 




Qc 


6 A 1 


Spreading factor for control symbols. 


45 




G /?' 


Spreading factor for i* code of A* user. 




r 


g^'XlA/c+u+if-D-i) R ece i ve vector associated with a time-slot. 






«*' 


Asynchronous delay (samples) for A* h user at the 








/71 th antenna. 


SO 


w 


6*' 


Channel impulse response length (chips). 




x {k} 


6 


Vector of symbols transmitted via i* code of A th user. 




*<*> 


6 jT 


Vector of symbols transmitted by A** 1 user. 


55 


X 




Vector of symbols transmitted by all users in a time-slot. 




X 


• A* 


Permuted vector of symbols transmitted by all users in a 








time-slot. 
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Claims 

1 . A method for detecting user symbols in a multi-rate multi-user environment, of a WCDMA communication system, 
the method comprising altering the order in which information symbols are estimated so as to group user symbols 

5 in a block-banded structure. 

2. A method as claimed in Claim 1 , wherein each user transmits control symbols and at least one group of data 
symbols, and wherein a convolution matrix < m )A is formed which is block-banded, each block including a number 
of columns relating to the number of groups of data or symbols transmitted by the users, and grouped according 

10 to the respective control symbols and in an order related to the order at which symbols arrive at a receiver. 

3. A method as claimed in Claim 2, wherein each block contains a control symbol for each user and one or more 
data symbols for each user. 

15 4. A method as claimed in Claim 3, wherein each block for n users comprises, sequentially, a control symbol for a 
first user i, one or more data symbols for the first user i, and, in turn a control symbol for each (i + 1) th user, and 
one or more data symbols for each (i + user. 

5. A method as claimed in any preceding claim, wherein the number of data symbols for each user within a block is 
20 determined by the number of groups of data transmitted by that user and the relative spreading factors of the 

control symbols and data symbols transmitted by that user. 

6. A method as claimed in Claim 5, wherein each user transmits synchronously a number N of control and data 
symbols, comprising one control symbol and N-1 data symbols, in an epoch, the method comprising shuffling the 

25 order to form a number N c of blocks, where N c is equal to the number of control symbols transmitted in a time slot, 

each block being of K c columns, and shuffling the blocks so that the order of the columns in each block is related 
to the time at which symbols arrive at a station. 

7. A method as claimed in Claim 1 , wherein each user transmits a plurality of groups of symbols wherein a convolution 
30 matrix ( m >A is formed which is block-banded, each block including a number of columns relating to the number of 

groups of symbols transmitted by the users, and grouped in an order related to the order at which symbols arrive 
at a receiver. 

8. A method as claimed in Claim 7, wherein each block for n users comprises, sequentially, a plurality of data symbols 
35 for a first user i, related to the number of groups of data symbols transmitted by that user, and one or more data 

symbols for each subsequent (i + 1 )* user, related to the number of groups of data symbols transmitted by each user 

9. A method as claimed in any preceding claim, wherein the number of data symbols for each user within a block is 
determined by the number of groups of data transmitted by that user and the relative spreading factors of the data 

40 symbols transmitted by that user. 

10. A method as claimed in Claim 9, wherein each user transmits synchronously a number N of data symbols, in an 
epoch, the method comprising shuffling the order selected to form a number N c of blocks, where N c is equal to the 
number of data symbols of a selected group transmitted in a time slot, each block being of Kc columns, and shuffling 

45 the blocks so that the order of the columns in each block is related to the time at which symbols arrive at a station. 

11. A method as claimed in Claim 10, wherein the selected group is the one having symbols of longer duration. 

12. A method as claimed in any preceding claim, substantially as hereinbefore described with reference to, and as 
so illustrated by, Figure 5 or Figure 6 of the accompanying drawings. 
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